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We report on charged hadron production in deuteron-gold reactions at ^/snn = 200 GeV. Our 
measurements in the deuteron-direction cover 1.4 < r\ < 2.2, referred to as forward rapidity, and in 
the gold-direction —2.0 < r/ < —1.4, referred to as backward rapidity, and a transverse momentum 
range pr = 0.5 — 4.0 GeV/c. We compare the relative yields for different deuteron-gold collision 
centrality classes. We observe a suppression relative to binary collision scaling at forward rapidity, 
sensitive to low momentum fraction (x) partons in the gold nucleus, and an enhancement at backward 
rapidity, sensitive to high momentum fraction partons in the gold nucleus. 

PACS numbers: 25.75.Dw 



Deep inelastic scattering of leptons on the proton re- vealed the proton's substructure of point-like parton con- 
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stituents This substructure, usually described quan- 
titatively as Parton Distribution Functions, evolves as 
one probes the proton at shorter wavelength or equiv- 
alently higher momentum transfer, Q 2 . Using the mea- 
sured quark and antiquark distribution functions and the 
DGLAP and BFKL [j| evolution equations, a strong 
increase in the gluon density is expected at high Q 2 
and small x (fraction of the proton momentum carried 
by the parton). Such an increase is indeed observed at 
HERA 4], suggesting that at sufficiently small x, gluons 
should overlap in space and time. This overlap should 
result in gluon fusion, and thus reduce the gluon density 
at low x and enhance it at larger x. This gluon fusion 
limits the achievable gluon density, leading to gluon sat- 
uration. This saturation is sometimes described as the 
formation of a Color Glass Condensate (CGC) 0. Gluon 
saturation is expected to be a larger effect in nuclei where 
the partons from different nucleons overlap as well. Sup- 
pression of low x partons in nuclei relative to nucleons 
has been experimentally observed and is referred to as 
nuclear shadowing However, this shadowing is often 
described in terms of modification of the leading-twist 
parton densities in nuclei 0. 

In 2003 the Relativistic Heavy Ion Collider (RHIC) 
collided deuteron and gold nuclei at ^snn = 200 GeV. 
At this energy, most hadrons with px > 2.0GeV/c arise 
from parton-parton interactions and can be used as a 
probe of nuclear partonic structure. Hadrons with px > 
2.0GeV/c at forward rapidity 1.4 < rj < 2.2 are sensitive 
to low x partons in the gold nucleus 0.001 < x < 0.03. 
Hadrons at backward rapidity —2.0 < rj < —1.4 are sensi- 
tive to high x partons in the gold nucleus 0.04 < x < 0.5. 
It has been predicted that gluon saturation at small x 
will suppress hadronic yields at forward rapidity [8( with 
the transverse momentum scale for the onset of the gluon 
saturation set by Q 2 [GeV 2 ] = 0.13 N coU e Xy @ for d + Au 
collisions at RHIC. Here A ~ 0.3 is determined from 
HERA data |l(J and N co u is the number of nucleon- 
nucleon inelastic collisions. Thus, for central collisions 
and within our forward rapidity coverage Q 2 is expected 
to be of order 2 — 4 GeV 2 and may have observable con- 
sequences. Novel hadron production mechanisms, such 
as quark recombination 11], can also impact the distri- 
bution of particles in the forward rapidity region. 

Results on charged hadron yields at forward rapidity 
from the BRAHMS experiment have shown a suppression 
of the yield of hadrons in central, compared to periph- 
eral, d+Au collisions [l^. At mid-rapidity, PHENIX has 
reported a modest enhancement of the yield of hadrons 
withp T > 1.5GeV/c 13]. This enhancement, generally 
referred to as the "Cronin effect" is often ascribed to ini- 
tial state scattering of the parton traversing the nucleus 
prior to the high Q 2 scattering ^4|. At backward ra- 
pidities (large x) , anti-shadowing and other effects of the 
surrounding nuclear medium (e.g. the EMC effect) 0] 
may compete, making predictions challenging. 



It is important to note that in the transverse mo- 
mentum range of this measurement, 0.5 < px < 
4.0GeV/c, hadron production is also sensitive to soft 
physics phenomena which are determined by coherent 
hadron-hadron interactions. In p + A reactions at lower 
energies soft hadron production shifts from forward to 
backward rapidity, with a larger shift for larger nuclear 
targets. Thus, at low px one may observe an increase 
(decrease) in hadron yields at backward (forward) rapid- 
ity which is not necessarily a reflection of changes at the 
partonic structure level. 

In this Letter, we present results from the PHENIX ex- 
periment [lfij ] on the ratio of hadron yields at forward and 
backward pseudorapidity for different centrality classes 
of d + Au collisions. PHENIX has two spectrometers 
designed for measuring muon production over the pseu- 
dorapidity range —2.2 < rj < —1.2 (backward spectrom- 
eter) and 1.2 < r] < 2.4 (forward spectrometer) yj|. The 
spectrometers start with a thick hadron absorber com- 
prised of 19 cm of brass and 60 cm of low-carbon steel 
between the collision point and active detectors along 
the beam axis, primarily to reduce hadronic background 
for muon measurements [TtI ]. After this material, the 
Muon Tracker (MuTr) detector, consisting of three sta- 
tions of cathode strip chambers, tracks charged particles 
in a magnetic field. The momentum resolution is 5% (for 
typical momenta in this analysis) and the absolute scale 
is known to better than 1%. Following the muon magnet 
backplate (30 (20) cm of steel in the forward (backward) 
spectrometer) there is a Muon Identifier (MuID) detec- 
tor. The MuID consists of five layers of planar drift tubes 
interleaved with layers of steel for further hadron absorp- 
tion (10 cm thick in the first two layers and 20 cm thick 
for the remaining layers). The layers are numbered 0-4, 
with 4 being the most downstream. The MuID is used 
to separate muons from hadrons and provide triggering 
capabilities. 

Although these spectrometers were designed to de- 
tect muons, they can also be used to measure charged 
hadrons via two independent methods. The first method 
is via the identification of hadrons which penetrate part 
way through the MuID, referred to as "punch-through 
hadrons." The second method is via muons from light 
mesons 7r, K which decay before interacting in the ab- 
sorber material. By measuring these decay muons, we 
can reconstruct the yield of their parent light mesons. In 
both of these methods the absolute yield of hadrons is 
difficult to determine due to uncertainties in the punch- 
through and decay probabilities. However, this small 
probability is independent of d + Au collision centrality, 
and thus not knowing the absolute yields does not affect 
the precision of measured ratios of hadron yields between 
the different classes of events. 

The punch-through hadron identification is achieved 
by studying particles that stop somewhere within the 
MuID before the last layer. For muons penetrating up 
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to layer 2 and 3 we expect the average momenta mea- 
sured in the MuTr of p = l.OGeV/c and p = 1.2GeV/c, 
respectively, corresponding to the average ionization en- 
ergy loss in traversing the spectrometer material. The 
reconstructed momentum distributions for particles stop- 
ping in layers 2 and 3 of the MuID are shown in Figure^ 
In addition to the expected muon peaks, there is a broad 
distribution extending to higher momentum which is the 
result of punch- through hadrons. These hadrons also suf- 
fer ionization energy loss up to the relevant layer and 
then suffer an inelastic collision in the MuID steel and 
do not penetrate further. We thus select a clean sam- 
ple of hadrons by demanding that a track stop in MuID 
layer 2 or 3 and have momentum more than 3<r away from 
the muon peaks. Muon contamination in our sample is 
estimated from simulations to be less than 5%. 
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FIG. 1: (color online), (left) The total momentum ptot 
measured in the MuTr without energy loss correction of all 
charged tracks penetrating to MuID layers 2 (grey) and 3 
(black), (right) Collision vertex distribution for events with 
muons at forward rapidity, corrected for the minimum bias 
collision vertex distribution. 

Another source of background for the punch-through 
hadrons is secondary particles produced from hadronic 
showering in the absorber. This background is reduced 
by requiring that the track point back to the primary 
collision vertex, as determined from the Beam-Beam 
Counter (BBC). The background from secondary par- 
ticles varies as a function of pr and is typically ~ I — 5% 
of the signal based on simulations. We also apply accep- 
tance cuts —2.0 < T) < —1.4 and 1.4 < -q < 2.2 in order 
to reduce background at small angles. 

Some hadrons will decay into muons before the ab- 
sorber, and the decay muons are then measured by 
the muon spectrometers. Muons can result from many 
sources including decays of n, K, D mesons, and J/ip- 
These particles have a finite decay probability Pdecay be- 
fore they reach the absorber 



Pdecay [P, L) — 1 — 



(1) 



where L ~ 41cm, is the distance from the collision vertex 



to the absorber; p, m and r are the momentum, mass and 
proper lifetime of the parent particle. 

Thus, collisions that occur far from the absorber will 
be more likely to produce muons from light meson de- 
cays than those that occur close to the absorber. Charm 
hadrons, however, due to their very short proper decay 

Lm 

lengths, e T -p < 1, will have minimal collision vertex 
dependence . The right panel of Figure^shows the colli- 
sion vertex distribution from events in which muons are 
detected at forward rapidity, corrected for the minimum 
bias collision vertex distribution. The large vertex de- 
pendence indicates a significant fraction of the muons are 
from pion and kaon decay. Using this distribution, we can 
separate the muons from pion and kaon decay from other 
contributions. The acceptance and efficiency vary by less 
than 5% over the z vertex range, which establishes an 
upper bound on the systematic error attributable to the 
subtraction of these non-signal contributions. It should 
be noted that the measured muon px is approximately 
15% lower on average than the parent hadron px, which 
is not corrected for in this analysis. 

The data set for this analysis was collected under two 
different trigger conditions. We recorded 67 x 10 6 min- 
imum bias triggers which required at least one hit in 
both the PHENIX forward 3.0 < r/ < 3.9 and backward 
—3.9 < rj < —3.0 Beam-Beam Counters (BBC) and a re- 
constructed vertex position within \z\ < 30 cm along the 
beam axis. The minimum bias trigger accepts 88 ± 4% of 
all inelastic d + Au collisions The second data set, 
sampling 5.3 x 10 9 minimum bias events, was collected 
with the MuID trigger which requires at least one track 
penetrating the first four layers of the MuID. 

We divide these events into four centrality classes 
based on the number of particle hits in the backward 
BBC counter covering —3.9 < rj < —3.0. Using a Glauber 
model 0| and simulation of the BBC, we determine the 
average number of binary collisions in each centrality 
class. The classes are categorized as follows: 60 — 88% 
((N coll ) = 3.1 ± 0.3), 40 - 60% ((N coll ) = 7.0 ± 0.6), 
20 - 40% ((N cM ) = 10.6 ± 0.7), and - 20% ((N coll ) = 
15.4 ±1.0). 

There is a correlation between having a particular 
physics process (for example the production of a high 
Pt hadron) and the BBC response. The BBC coverage 
in pseudorapidity is well separated from the muon spec- 
trometers so the correlation is not predominantly due to 
jet fragmentation, but rather an underlying event correla- 
tion. We have studied this effect in detail using proton- 
proton and d + Au data, and in simulations, and have 
accounted for this correlation bias. The bias correction 
factors we apply range from 0-7% depending on the cen- 
trality category and the physics process. The systematic 
errors on these corrections are less than 4%. 

The nuclear modification factor R cp is defined as the 
ratio of the particle yield in central collisions to the par- 
ticle yield in peripheral collisions, each normalized by 
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FIG. 2: (color online). R cp as a function of pr at forward 
rapidity (squares) and backward rapidity (circles) for different 
centrality classes. 
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FIG. 3: (color online). R cp as a function of r\ for 1.5 < pr < 
4.0GeV/c for different centrality classes. 



point systematic errors result from sensitivities to analy- 
sis cuts and are 5 — 10%. They are added in quadrature 
with the statistical errors and shown as error bars. 

It is notable that our two measurement methods have 
different sensitivity to different hadrons. The particle 
composition (n/K/p ratio) of the observed sample is 
modified relative to the particle composition at the colli- 
sion vertex due to species-dependent nuclear interaction 
cross-sections affecting the punch-through hadrons and 
due to species-dependent decay lifetimes affecting the 
hadron decay muons. Both effects enhance the kaon con- 
tribution to our R cp measurements. The uncertainty on 
our charged hadron R cp values introduced by this effect 
is estimated to be less than 4% by calculating the differ- 
ence between the kaon R cp and inclusive char ged particle 
R cp determined by PHENIX at mid-rapidity 18] . 

We observe that R cp shows a suppression at forward 
rapidity that is largest for the most central events. The 
opposite trend is observed at backward rapidity where 
R cp shows an enhancement that is also largest for the 
most central events. We observe a weak pr dependence 
with slightly smaller R cp values at lower px- We observe 
a clear pseudorapidity dependence at forward rapidity 
with R cp dropping further at larger rj values. Within 
our current uncertainties we are unable to discern any 
pseudorapidity dependence at backward rapidity. 

In Figure 4 we compare results from the BRAHMS 
experiment |l2l| with our results at forward rapidity. The 
PHENIX data and the BRAHMS data are in agreement 
within systematic uncertainties. 



the average number of nucleon-nucleon binary collisions 
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The hadron R cp , using the most peripheral centrality 
class (60 — 88%) for normalization, is shown in Figure [21 
as a function of px at forward and backward rapidities. 
The results from the punch-through hadron (PTH) and 
hadron decay muon (HDM) techniques are both shown 
and are in quite good agreement. We also show the re- 
sults integrated over 1.5 < pr < 4.0GeV/c as a function 
of pseudorapidity in Figure [3J 

There are two types of systematic uncertainties in 
our analysis. Common systematic errors which move 

all data points up and down together include the error 

i 

— (10.8% for the most central bin), the cen- 



N p ' 



trality bias correction factors (4%), and the centrality- 
dependent tracking efficiency (4%) determined by em- 
bedding Monte Carlo particles in real data. Common 
systematic errors are shown as a black bar. Point-to- 
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* BRAHMS ii=2.2 
BRAHMS ri=3.2 
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FIG. 4: (color online). PHENIX R cp as a function of pr at 
forward rapidities shown as the average of the two methods. 
Note that the BRAHMS results are for negative hadrons at 
77 = 2.2,3.2 and their centrality ranges (0 - 20%/60 - 80% 
and 30 — 50%/60 — 80%) are somewhat different from ours. 



The suppression of hadron yields relative to binary col- 
lision scaling at forward rapidity is expected from ini- 
tial state nuclear effects. However, detailed comparisons 
with various theoretical approaches is necessary in order 
to discriminate between different models. In particular 
the lack of a strong pt depedence at both forward and 
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backward rapidities must be understood as the physics 
processes transition from "soft" to "hard" physics scales. 

To summarize, we observe a suppression in hadron 
yields relative to binary collision scaling at forward ra- 
pidities and an enhancement at backward rapidity for 
central relative to peripheral d+Au reactions at ^snn — 
200 GeV. The forward rapidity suppression is in qualita- 
tive agreement with the expectation of shadowing and 
saturation effects in the small x region in the gold nu- 
cleus. However, other physics effects must also be con- 
sidered in understanding the full px and r\ dependence. 
The source of the backward rapidity enhancement, and 
the possible contribution of anti-shadowing of large x par- 
tons, has yet to be understood. 
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